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1 Introduction 

In this contribution, we report on the dynamical coupled- 
channel analysis being pursued at the Excited Baryon 
Analysis Center (EBAC) of Jefferson Laboratory. EBAC 
was established in January, 2006. Its objective is to extract 
the parameters associated with the excited states (N*) of 
the nucleon from the world data of meson production re- 
actions, and to also develop theoretical interpretations of 
the extracted N* parameters. 

Since A^* states are unstable, their structure must cou- 
ple with the reaction channels in the meson production re- 
actions. To determine correctly the spectrum of N * states, 
an analysis of the meson production data must account 
for the coupled-channel unitary condition. The extracted 
N* parameters can be interpreted correctly only when the 
reaction mechanisms in the short-range region where we 
want to map out the N* structure have been accounted 
for. To meet these two crucial requirements, the dynam- 
ical coupled-channel model developed in Ref.p] is being 
applied at EBAC to analyze the data of n, tttt, r}, K, and 
ui production. It involves sloving the following coupled in- 
tegral equations 

T a ,f3(p a ,Pf3\E) = Va,p{p a ,pfi) 
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the near future, one hopes to relate -TV*,/? to Lattice QCD 
(LQCD) calculations. 

If we take the on-shell approximation, Eq.{T]) is re- 
duced into an algebraic form of K-matrix models 0,5. 6, 
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where p(p 7 ) is an appropriate phase space factor. Quali- 
tatively speaking, the use of the on-shell model based on 

piPaiPp, E) of Eq.([3]) is to avoid an explicit treatment 
of the reaction mechanisms in short range region where 
we want to map out the quark-gluon sub-structure of A^* 
states. Thus the N* parameters extracted by using Eq.(p} 
can be more directly interpreted in terms of the quark- 
gluon sub-structure of N* . 

In section 2, wc briefly review the dynamical coupled- 
channel model of Ref.p]. The recent results from EBAC 
are reported in section 3. Section 4 gives an outlook of 
EBAC. 



2 Dynamical Coupled-Channel Model 

Within the formulation of Ref.[T], Eqs.lQ])-© are derived 
from the following Hamiltonian 



E — M* 



(2) 



H = Hn 



Hr 



where a, /?, 7 = 7JV, ttN, ?yA^,and 7T7rA^ which has irA, pN, aN 
resonant components, G 1 (p, E) is the propagator of chan- 
nel 7, v a> j3 is defined by meson-exchange mechanisms, and 
En*,p is related to the quark-gluon sub-structure of N*. 
At the present time, it is reasonable to interpret -T/v-,/3 in 
terms of hadron structure calculations with effective de- 
grees of freedom, such as the constituent quark model|2] 
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The non-resonant interactions ^22,^23,^33 are derived 



and the model[3] based on Dyson-Schwingcr Equations. In from phenomenological Lagrangians by using an unitary 
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Fig. 1. Basic mechanisms of the Model Hamiltonian defined 
in Eqs.(4)-(6). 
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Fig. 2. Meson-baryon (MB) interaction mechanisms of «2,2 of 
Eq.(6). 
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Fig. 4. One-particle-exchange interactions Z^ 7vA (E), 
KnIa ^ ZfJ^ ofEq.(9). 
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Fig. 5. Examples of mechanisms included in Z$ B M , B , of 
Eq.(9) 



(d) 



(e) 



The resonant amplitude in Eq.0 is defined by 

ta, b {E)= £ rl : JE)lg(E)}^r N; , b (E), (10) 
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Fig. 3. Examples of non-resonant mechanisms of vmn,wwN 
with M — n or 7 (denoted by long-dashed lines). Mi denotes 
the intermediate mesons (ir,p,ui). 

transformation method [9lll0j. The interaction V22 is de- 
fined by the tree diagrams, as illustrated in Figj2] Exam- 
ples of non-resonant mechanisms of of V23 are illustrated in 
Fig[3J The more complex W33 is not included in the current 
analysis at EBAC. 

By using the projection operator techniques, one can 
cast Eqs.([T])-([2]) for 2^2 meson-baryon (MB) transition 
amplitudes into the following form 



where the dressed vertex is 

r N *, a (E) = r N . ;a + J2rN*, b G b (E)t bia (E), (11) 
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and the N* propagator is defined by 

U,J (E)=E- M?5u - £ f N * , a G a {E)r%, >a . (12) 



T a , b (E) = t a , b (E)+t*JE) 
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where a,b,c= jN, irN, rjN, irA, pN, aN . The non-resonant 
amplitude in Eq.([7]) is defined by 

t a ,b(E) = V a , b + ]T V a , c G c (E)t c , b (E) , (8) 

C 

where the driving term is 

V a , b (E) = va.b + Z { a E b \E) + Z§{E). (9) 

The effects due to the three-body unitarity cuts arc in- 
cluded in Z ( a E b \E) and Z^E), as illustrated in FigsgHSl 



Here we emphasize that the second term of the dressed 
vertex Eq. (|ll[) is an necessary consequence of the unitar- 
ity condition. This term describes the meson cloud effects 
on the N* excitations. See Rcf. [8] for a more detailed dis- 
cussion on this point. 

To solve Eq.©, we need to handle the singular struc- 
ture of Z^ and Z^ b . Their matrix elements diverge within 
the moon-shapes region in FigJB] and thus have singular 
structure illustrated in Fig [7] We overcome this difficulty 
by using the Spline- function method developed in Refs. fTTl 

Eg. 



3 recent Results 

The computer code for solving the dynamical coupled- 
channel equations given in section 2 has been developed. 
The numerical details, in particular on the use of Spline- 
function methods, are explained in Ref.pQ. 
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Fig. 6. Logarithmically divergent moon-shape regions of 
the matrix elements of n ^{p' ,P, E) (solid curves) and 
z( p n]iia(p'>P> E ) (dashed curves). 
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Fig. 7. Matrix elements of the one-particle-exchange term 
zf2 wA {k,k',E) for L = L' = 1, J = 5/2, T — 1/2 at k' = 300 
MeV/c and E = 1.88 GeV. 



Our first task is to determine the parameters of the 
hadronic interactions of Hamiltonian defined by Eq.®. 
This has been done[T5] by fitting the irN scattering data. 
One or two bare N* states in each of S, P, D, and F 
partial waves are included to generate the resonant am- 
plitudes, defined by Eq. (fT0|) . in the fits. The parameters 
of the model are first determined by fitting as much as 
possible the empirical ttN elastic scattering amplitudes 
of SAID [4] up to 2 GeV. Some of our fits are shown in 
Figs BH 




Fig. 8. Fit to the I = § Re(T^ N>nN ) of SAID [4]. 



We then refine and confirm the resulting parameters 
by directly comparing the predicted differential cross sec- 
tion and target polarization asymmetry with the origi- 
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Fig. 9. Fit to the I = \ Im(T nNt7rN ) of SAID [4], 
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Fig. 10. Differential cross section for several different center 
of mass energies. Solid curves correspond to our model while 
blue dashed lines correspond to the SP06 solution of SAID [4]. 
All data have been obtained through the SAID online applica- 
tions. Ref. [4]. 



nal data of the elastic 7r ± p — > 7r ± p and charge-exchange 
tt~p — > 7r°n processes. Typical results are shown in FigfTUl 
The predicted total cross sections of ttN reactions are also 
in good agreement with the data, as shown in the left hand 
side of FigQl] Our predictions of the partial total cross 
sections to each channel, as given in the right hand side of 
FigfTTl need to be refined by also fitting the irN — > ititN 
data. A combined fit to both the data of ttN elastic scat- 
tering and 7rA — ► ttttN is now being pursued at EBAC. 
FigH21 is a result from this very challenging task, show- 
ing the importance of the 7T7tA cut in determining the 
predicted invariant mass distribution of ttN — > ttttN reac- 
tion. Similar pronounced effects on 7./V — > 7T7tA reactions 
have been presented in Ref. [T]. 
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Fig. 11. Predicted -n p total cross sections 

With the hadronic parameters determined, we arc now 
moving to analyzing the data of electromagnetic produc- 
tion of 7t and 7T7T. Here the only parameters to be de- 
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Fig. 14. Coupled-channel effects on the ~fN — > tyN in the A 
region 



Fig. 12. The -kN invariant mass distribution of nN 
reaction at W = 1880 MeV. 
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termined are the bare jN — > N* vertex interactions. In 
Fig |131 we show the preliminary results for — ► Tr + n. In 
Fig. 14 we show the coupled-channel effects on the 7^V — ► 
irN reactions in the A excitation region. More detailed 
results will be presented in Ref.|14|. 



4 Outlook 

The dynamical coupled-channel model developed in Ref. [T] 
is being applied at EBAC to analyze the world data of 
meson production reactions. The analysis of ir and irir 
production data is proceeding well with some encourag- 
ing preliminary results. We have also started to analyze 
the data of rj, K , and uj production. Here the main task is 
to extend the model to include the interactions associated 
with the KA, KS, and uiN channels. 

In parallel, we are also developing an analytical contin- 
uation method to extract the resonance poles and residues 
from the partial wave amplitudes predicted by the con- 
structed dynamical coupled-channel model. Furthermore, 
we are investigating the extent to which the resonance 
poles and bare N* parameters extracted from our analy- 
sis can be related to the current hadron structure calcu- 
lations; in particular the Lattice QCD calculations. 
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Fig. 13. Fit to the 7p -> n + p up to W = 2 GeV 



